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1.  INTRODUCTION 


This  chapter  discusses  the  basic  issues  of  the  thesis.  We  begin  by  describing  fracture 
mechanics  and  discussing  why  it  is  important.  We  continue  with  a  brief  explanation  of  linear  elastic 
fracture  mechanics  and  go  on  to  discuss  scaling  and  why  it  is  important.  Then  we  outline  the  size 
predictions  that  are  obtainable  from  LEFM  for  three  particular  situations.  After  this  we  define  the 
objectives  of  the  present  work  and  finally  we  give  a  section-by-section  listing  of  the  topics  to  be 
covered  herein. 


One  way  to  discuss  fracture  mechanics  is  to  realize  that  objects  are  made  from  materials 

that  are  not  perfect.  These  imperfections  or  defects  can  cause  many  things  to  happen.  The  thing 

that  is  of  importance  here  is  that  these  defects  can  cause  a  local  stress  increase  in  the  object. 

The  amount  of  this  stress  increase  is  related  to  the  geometry  of  the  defect  itself  and  the  object 

within  which  it  is  contained.  As  the  geometry  of  the  defect  becomes  more  acuminated,  the  local 

stress  begins  to  approach  infinity,  i.e.  there  is  a  stress  singularity.  This  analytical  result  creates  a 

conflict  with  experience  since  there  is  common  knowledge  of  items  containing  sharp,  crack-like 

features,  which  indeed  do  not  break  when  exposed  to  the  slightest  of  loadings.  Fracture 

mechanics  addresses  itself  to  this  inconsistency.  A  more  formalized  definition,  based  upon 

material  defects,  is  given  by  Kanninen  and  Popelar  [1],  p.  89. 

"Fracture  mechanics  is  an  engineering  discipline  that  quantifies  the  conditions 
under  which  a  load-bearing  solid  body  can  fail  due  to  the  enlargement  of  a  dominant 
crack  contained  in  that  body." 


Some  indication  as  to  why  one  might  be  interested  in  fracture  mechanics  is  highlighted  in 
reports  by  Duga  et  al.  [2]  and  [3],  These  references  cover  a  study  conducted  by  Battel 
Laboratories  for  NBS  in  which  the  cost  of  fracture  in  the  U.S.  is  placed  at  1 19  billion  dollars  (1982 
$)  per  year  which  is  approximately  4  per  cent  of  our  Gross  National  Product.  This  study  also 
indicated  that  80  per  cent  cf  that  total  is  associated  with  efforts  to  prevent  fracture,  as  opposed  to 
simply  replacing  broken  parts.  The  study  states  that  if  the  newest  fracture  mechanics 


technologies  were  applied  to  manufacturing  techniques,  then  $35  billion  of  this  total  could  be 
saved.  Based  on  these  numbers  one  can  see  that  it  would  be  economically  prudent  to  be  able  to 
better  understand  fracture  mechanics. 


The  most  fundamental  and  oldest  portion  of  fracture  mechanics  is  that  dealing  with  the 
response  of  a  linear  elastic,  or  brittle  material,  linear  elastic  fracture  mechanics,  LEFM.  An 
enjoyable  history  of  its  development  is  given  by  Kanninen  and  Popelar  [1],  LEFM  is  based  upon 
the  use  of  a  parameter  called  a  stress  intensity  factor,  K,  to  predict  fracture.  The  stress  intensity 
factor,  K,  is  defined  to  be  a  coefficient  of  the  singular  stress  field.  In  general,  the  form  of  K,  as 
shown  by  Broek  [5]  and  others,  is, 

K=oV^Y(^.  (1.1) 


where  G  is  the  nominal  applied  stress,  or  the  far-field  stress.  The  a  is  the  size  of  the  defect  or  half 
crack  length.  And  V(^)is  a  geometry  function  describing  the  shape  of  the  part;  it  is  based  upon 

the  ratio  of  defect  size  to  the  nominal  specimen  size  or  width,  W.  LEFM  says  that  when  this  K 
value  reaches  a  particular  level,  the  material  fractures.  The  result  is  that  by  computing  the  stress- 
intensity  factor  for  a  particular  load,  and  object  geometry,  one  can  determine  if  a  fracture  could 
occur  by  comparing  the  K  value  to  a  minimum  critical  stress-intensity  factor  «ic  for  the  material! 
However,  this  requires  an  elastic  material  response,  that  is  to  say,  there  should  be  no  plastic 
deformation  in  the  material.  In  applications,  however,  LEFM  is  used  if  any  region  of  plastic 
deformation  is  confined  to  some  relatively  small  distance  around  the  defect. 


t  The  parameter  K  is  defined  for  three  modes  of  loading,  as  shown  in  Rolfe  and  Barsom 
[4],  «icis  the  minimum  critical  stress  intensity  to  occur  for  mode  one  or  the  crack  opening  mode, 
in  plane-strain. 
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As  with  most  technologies,  the  usefulness  of  LEFM  is  greatly  enhanced  if  it  can  be 
applied  to  a  broad  range  of  specimen  sizes.  The  ability  to  change  with  size  or  scale  and  still 
maintain  its  applicability  is  a  basic  requirement  in  predictive  engineering  methods.  This  ability  to 
withstand  scaling  allows  a  technology  to  be  applied  from  situation  to  situation,  and  in  particular  it 
allows  a  technology  to  be  applied  from  a  small  scale  laboratory  testing  situation  to  an  application 
within  the  real  world.  By  way  of  illustration,  consider  the  difficulties  of  doing  a  full-scale  test  on  the 
wing  of  a  Boeing  747. 

This  study  will  look  at  the  question  of  scaling  in  LEFM,  and  in  particular,  the  effects  of 
scaling  on  two  distinct  types  of  loading  situations.  The  first  is  monotonic  loading,  basically  used  to 
determine  the  maximum  load  that  a  particular  geometry  can  withstand  or  the  critical  stress  intensity 
factor  for  a  material.  And  the  second  situation  is  that  where  the  loading  oscillates,  or  cycles.  The 
cyclic  type  of  loading,  usually  called  fatigue,  is  important  because  it  determines  the  useful  lifetime 
of  a  particular  part  at  a  load,  or  conversely,  the  load  able  to  be  endured  for  a  required  time  span. 

The  monotonic  scaling  question  can  be  broken  out  into  changes  in  two  separate 
dimensions,  one  being  how  wide  the  part  is  or  its  planar  size,  which  we  will  call  in-plane  scaling. 
This  is  represented  by  Figure  t  .1 .  As  can  be  seen  in  this  figure  both  width  and  defect  size  scale. 
In  a  similar  fashion  in  Figure  1 .2  is  shown  scaling  in  the  remaining  dimension,  or  out-of-plane 
scaling.  Note  that  in  this  figure,  the  frontal  view  of  the  specimen  does  not  change,  but  its 
thickness  does.  On  the  other  hand  for  scaling  under  cyclic  loading,  we  consider  scaling  without 
in-plane  or  out-of-plane  separations.  Therefore,  the  scaling  here  might  be  thought  of  as  an 
overall-scaled  or  total-scaled  situation.  This  geometry  is  represented  in  Figure  1.3. 


Now  we  go  on  to  discuss  what  LEFM  predicts  as  the  scaling  effects  in  each  of  these  three 
situations.  Addressing  first  the  in-plane  scaling  effect  for  the  monotonic  loading  situation,  if  one 
writes  (1.1)  for  each  of  the  two  sizes  of  Figure  1.1.  one  gets, 

Ki*-Oi*V«wY(§)  and  K2*  =  02*Vx^Y(^.  (1.2) 

Here,  the  superscript  *  indicates  a  critical  or  fracture  value.  LEFM  assumes  that  in  both  sizes 
failure  occurs  at  the  same  critical  K  value  if  the  specimens  are  made  of  the  same  material.  This 
allows  one  to  equate  Ki*  and  K2*  in  (1.2)  and  obtain  the  expression  for  the  strength  ratio  as, 

-21 (1.3) 
ct2 

What  (1.3)  says  is  that  as  the  cracked  specimen  gets  bigger  (X>1),  the  stress  at  which  it  fails,  or  its 
strength,  gets  smaller  relative  to  its  strength  for  its  initial  size. 

We  now  turn  our  attention  to  the  LEFM  prediction  for  out-of-plane  scaling  in  the 
monotonic  situation,  i.e.  that  demonstrated  by  Figure  1 .2.  Doing  a  similar  operation  as  above  for 
out-of-plane  scaling,  one  sees  that  there  is  no  explicit  statement  of  B  in  (1.1),  therefore,  there  is 
no  difference  between  the  critical  strength  of  the  thick  or  the  thin  situation.  This  is  the  accepted 
situation  for  perfectly  brittle  materials,  or  materials  with  limited  ductility.  However,  this  absence  of 
thickness  effect  changes  when  the  material  exhibits  a  degree  of  ductility.  First,  consider  a  thin 
specimen,  or  one  that  is  in  a  state  of  plane  stress,  that  is  where  there  exists  significant  stresses  in 
two  principal  directions  while  stresses  in  the  third  direction  are  negligible.  As  the  loading  of  this 
specimen  is  increased,  the  differences  between  these  principal  stresses  also  increases,  thereby 
increasing  the  shearing  stress.  Large  shearing  stresses  imply  yielding,  and  this  yielding  blunts 
the  crack  tip,  which  lowers  the  stress  intensification.  A  lower  stress  intensification  requires  more 
loading  to  cause  fracture,  therefore  the  fracture  strength  goes  up.  In  contrast,  consider  the 
situation  arising  when  a  thick  specimen  is  loaded.  Here,  there  exists  a  state  of  plane  strain,  i.e.  the 
stresses  of  the  three  princ:?al  directions  are  not  zero.  As  the  applied  load  is  increased,  the 
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differences  among  these  stresses  do  not  rise  as  quickly,  and  therefore  the  shearing  stress  is 
lower,  as  is  the  tendency  for  yielding.  Since  yielding  is  not  available  to  reduce  the  stress 
intensification,  the  thick  specimen  fractures  at  a  loading  less  than  for  the  thin  specimen.  In  light  of 
this  explanation,  the  critical  value  of  applied  load  is  expected  to  increase  with  reduced  thickness, 
hence  the  critical  value  of  the  stress  intensity  is  expected  to  increase  with  reduced  thickness. 


Finally,  we  turn  our  attention  to  the  cyclic  loading  situation.  The  equation  of  importance  in 
that  situation  is  known  as  the  Paris  Law, 


-jrr  =  C  (AK)m  for  AKmin  <AK<  AKmax  • 


This  empirical  equation,  based  on  a  data  reduction  scheme,  shows  that  the  crack  growth  per 
cycle,  da/dN,  is  a  function  of  the  change  in  the  stress  intensity  as  the  loading  cycles  from 
maximum  to  minimum,  AK,  and  two  material  parameters,  C  and  m.  What  this  expression  and  others 
like  it  attempt  to  do  is  predict  the  growth  rate  of  the  defect  from  an  initial  size  to  the  point  at  which 
the  crack  length,  a,  is  critical;  i.e.  a  gets  large  enough  so  that  K  of  (1 .1)  exceeds  the  critical  value. 
Equation  (1 .4)  has  a  useful  range  based  upon  the  change  in  stress  intensity.  That  is  to  say,  below 
some  "threshold"  AKmjn,  there  is  typically  no  crack  growth,  at  least  (1 .4)  does  not  apply  for 
whatever  growth  there  is,  and  above  AKmax  there  is  catastrophic  failure.  The  usual  means  of 
employing  (1 .4)  is  to  convert  it  to  an  expression  for  lifetime,  or  number  of  cycles,  N,  from  an  initial 
crack  size,  aj,  to  size  at  fracture,  af .  Rearranging  (1 .4),  substituting  an  expression  for  AK,  and 
preforming  the  necessary  integration  gives. 


where  AO  is  the  difference  of  extremes  of  the  applied  stress  during  the  cycle. 
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Essential  for  (1.5)  to  be  useful  in  practice  is  the  assumption  in  LEFM  that  C  and  m  are 
independent  of  size  or  scale.  Further,  for  successful  predictions  using  (1 .5),  there  is  really  an 
implicit  assumption  that  the  range  of  applicability  of  (1.4),  AKmjn  to  AKmax,  is  also  independent  of 
scale.  Under  these  conditions  it  is  interesting  to  examine  what  (1.5)  has  to  say  about  the  cyclic 
lives  of  small  and  large  specimens,  Ns  and  N|,  respectively.  Puttick  and  Atkins  (6)  discuss  this  and 
they  predict, 

*,nV2-,).  11.61 


Equation  (1.6)  assumes  af/W  and  af/W  are  the  same  in  both  sizes,  as  is  ACT.  Provided 
m>2,  the  life  of  the  small  specimen  is  predicted  to  be  greater  than  that  of  the  big. 


We  are  now  in  a  position  to  better  define  the  objectives  of  this  work.  The  goal  here  is  to 
attempt  to  do  a  critical  appraisal  of  how  well  these  accepted  predictions  from  LEFM  actually  work. 
We  feel  that  one  of  three  options  is  likely  to  occur.  The  first  is  that  LEFM  does  indeed  work  and, 
therefore,  is  a  viable  technology  and  needs  no  improvement.  The  second  possibility  is  that  the 
predictions  are  marginal,  that  is  to  say,  that  they  produce  the  right  general  trends  but  are  off  less 
than  an  order  of  magnitude.  This  would  imply  that  the  technology  is  not  without  redeeming 
features  but  needs  some  reworking  to  improve  its  capabilities.  The  last  possibility  is  that  LEFM,  as 
a  predictive  technology  simply  doesn't  work.  This  would  imply  that  it  is  not  even  correct  in  trends 
or  magnitudes  of  numbers.  To  handle  these  last  two  outcomes  is  beyond  the  scope  of  the 
present  work.  To  fully  reconcile  the  complete  body  of  data  obtained  over  the  last  thirty  to  forty 
years  under  these  circumstances  would  be  a  truly  Herculean  task,  though  some  directions  along 
which  one  might  proceed  in  this  event  can  be  suggested. 


The  remainder  of  the  thesis  is  organized  as  follows.  The  next  chapter  is  a  literature  review 
for  pertinent  data.  Following  this  is  a  description  of  a  series  of  monotonic  loading  experiments 
after  which  we  describe  a  series  of  cyclic  loading  experiments.  The  fifth  section  contains  some 
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2.  LITERATURE  REVIEW 


In  this  chapter  we  present  the  results  of  a  review  of  the  literature  pertinent  to  the 
questions  of  Chapter  1 .  We  begin  by  discussing  the  literature  that  addresses  in-plane  scaling, 
and  then  go  on  to  articles  on  thickness  scaling.  Finally,  we  conclude  by  looking  at  some  papers 
that  address  the  question  of  scaling  in  cyclic  loading  situations. 

With  respect  to  in-plane  scaling,  a  recently  completed  local  work,  Sinclair  and  Chambers 
[7],  comments  upon  how  well  the  physical  evidence  compares  to  predictions  based  upon  LEFM. 
After  an  extensive  review  of  the  open  literature,  the  authors  present  six  plots  of  normalized 
strength  versus  scaling  factor  which  summarize  their  results  for  plane-strain  or  plane-stress  and 
material  response  of  brittle,  brittle-ductile,  or  ductile  .  A  typical  curve  is  reproduced  in  Figure  2.1 . 
This  particular  figure  shows  the  plane-strain  brittle  response  for  a  number  of  reported  articles.  The 
LEFM  prediction  of  equation  (1.3)  is  shown  by  the  solid  line,  and  the  points  represent  reported 
results,  large  points  indicating  many  experiments.  As  can  be  seen  in  this  figure,  the  agreement 
between  the  prediction  and  the  experimental  results  is  not  good.  To  quantify  this  statement, 
somewhat,  LEFM  is  found  to  give  predictions  that  are  not  within  ±  10%  of  actual  results  over  80% 
of  the  time  when  the  scale  factor,  X,  is  greater  than  3,  and  1 00%  of  the  time  when  X  >  7.  The  article 
concludes  by  stating  that,  although  the  results  are  trendwise  correct,  they  tend  to  be  too  simplistic 
and  that  the  net  result  is  a  need  to  test  at  various  specimen  sizes. 

Now,  we  turn  our  attention  to  out-of-plane  scaling.  A  search  of  the  open  literature  was 
undertaken  for  results  of  experiments  scaling  thicknesses  of  otherwise  identical  specimens.  To 
be  compatible  with  the  above  article  [7],  many  of  the  same  guidelines  are  used  when  reporting  the 
results.  Comparison  is  made  only  within  a  single  source  and  only  to  the  extent  that  identical 
geometries  (center-cracked  tension, edge-notch  bending,  etc.), materials  and  environment  are 
used.  In-plane  dimensions  are  required  to  remain  constant  for  various 


thicknesses  used,  although  a  small  variation  of  crack  length,  a,  is  permitted;  i.e.,  A(a/VJ)  <  ±  1 0%. 
Crack  or  notch  acuity  must  be  sharp  as  determined  by  the  ratio  of  defect  radius  of  curvature,  p,  to 
a;  i.e.,  p/a  <  0.1 :  fatigue  precracking  is  preferred.  The  material  response  is  categorized  into  three 
regimes,  based  upon  the  relative  plastic  zone  size,  ry/a;  i.e. : 


ry/a  <  0.02 
0.02  <  ry/a  <  0.05 
0.05  <  ry/a 


brittle; 


brittle-ductile,  and; 


ductile. 


Relative  plastic  zone  size  is  computed  by  using: 


ry/a-—  f— 
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for  plane  stress,  or 


for  plane  strain,  or, 


between  the  two. 


Here,  Oy  is  the  yield  strength  of  the  material  at  similar  conditions.  The  classification  according  to 

material  response  is  done  since  a  brittle  response  affords  the  best  comparison  with  LEFM,  and 
the  brittle-ductile  response  is  the  natural  place  into  which  extension  of  the  LEFM  prediction  is 
logical.  Strength  is  taken  as  the  nominal  applied  stress  at  the  onset  of  Mode  I  crack  growth.  This 
can  be  inferred  from  energy  release  rate  (G),  stress  intensity  (K),  or  nominal  net  stress,  although 
the  preferred  information  is  critical  load  and  specimen  geometry.  If  the  only  information  given  is  at 
maximum  loading,  it  was  used,  otherwise  the  values  at  5%  offset  or  at  pop-in  were  preferred.  This 
is  done  to  get  information  at  the  onset  of  crack  growth,  when  plastic  rioformafion  shcu'd  be 
minimal. 


After  review  of  some  articles,  it  became  apparent  that  simply  reporting  strength  versus 
thickness  would  not  work.  By  way  of  illustration  consider  two  experiments:  in  one,  the  in-plane 
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dimension  is  one  inch  and  the  thickness  varies  from  1/8  inch  to  3/4  inch.  In  the  other,  the 
thickness  varies  in  the  same  range,  but  now  the  major  in-plane  dimension  is  12  inches.  The  point 
here  is  that  the  same  thickness  in  one  experiment  will  clearly  be  in  a  plane-strain  regime,  while,  in 
another  experiment,  that  thickness  is  in  the  plane-stress  regime.  To  overcome  this  problem,  we 
decided  to  characterize  results  based  upon  a  dimension  that  would  essentially  define  the  plane- 
stress  to  plane-strain  transition  for  all  test  specimens.  To  accomplish  this,  we  somewhat  arbitrarily 
use  the  ratio  of  thickness  to  crack  length,  specifically  2B/a  as  the  transition  parameter.  This  is 
comparable  to  the  ASTM  standard  E399  [8]  which  uses  the  ratio  of  thickness  to  width  to  define 
plane-strain  response.  If  the  ratio  2B/a  equals  one  for  a  configuration,  then  the  specimen  is  taken 
to  be  at  the  transition  from  the  plane-stress  and  plane-strain  loading  regions.  If  the  ratio  is  two  or 
better,  then  the  response  is  plane-strain,  while  if  the  ratio  is  less  than  0.2,  the  response  is  plane- 
stress.  The  reported  strength  is  normalized  by  the  strength  for  a  specimen  with  a  2B/a  ratio  equal 
to  one.  If  such  a  specimen  was  not  actually  reported,  then  the  data  is  interpolated  if  possible, 
extrapolated  where  not,  to  the  proper  ratio.  The  basis  for  the  interpolation/  extrapolation  is  a 
straight  line  on  a  log-log  plot.  Clearly,  this  is  not  the  only  way  to  handle  this  scaling  situation,  but, 
hopefully,  it  represents  one  rational  way  of  doing  it. 

In  the  course  of  this  phase  of  the  work,  in  excess  of  1 20  articles  were  reviewed  for  data. 

Of  these,  only  58  [9  to  67]  contained  enough  data  to  allow  comparison  of  relative  strength  to 
relative  thickness  .  The  remaining  62  articles  [68  to  130]  were  such  that,  although  they  address 
the  topic  of  out-of-plane  scaling,  they  did  not  contain  enough  information  to  either  calculate 
relative  strength  or  size,  or  to  compute  the  response  of  the  material.  Such  articles  are,  though, 
included  for  completeness.  Although  this  data  search  is  probably  not  all-inclusive,  it  should  give  a 
feeling  for  the  trend  of  the  reported  experimental  results.  The  detailed  information  that  could  be 
obtained  from  these  articles  is  contained  in  tables  in  Appendix  A.  This  same  information, 
classified  as  to  material  response,  is  plotted  in  Figures  2.2, 2.3,  and  2.4. 


Figure  2.2  presents  the  dependence  of  strength  upon  thickness  when  ry/a  <  0.02,  i.e. 
brittle  response,  which  should  most  closely  fit  the  standard  LEFM  prediction.  If  this  response 
determination  where  checked  again  using  ry/B  <  0.02,  it  is  possible  that  in  some  cases  the 
thinnest  specimens  could  no  longer  be  classified  as  brittle.  However,  in  order  to  include  more 
data,  we  restrict  ourselves  to  the  ry/a  criteria.  As  can  be  seen,  most  of  the  data  groups  about  a 
scaling  factor  of  one,  indicating  that  most  of  the  testing  is  done  in  the  plane-stress  transition 
region.  The  LEFM  prediction  for  this  response  is  a  horizontal  line  at  a  relative  strength  of  one.  As 
can  be  seen,  the  results  oscillate  above  and  below  this  prediction.  Relatively  few  experiments  are 
actually  represented  on  this  curve,  and  the  range  in  any  single  experimental  study  is  rather  small. 
Presented  for  completeness  are  similar  reported  data  for  material  responses  in  the  brittle-ductile 
and  ductile  regimes,  plotted  in  Figures  2.3  and  2.4,  respectively.  Note  in  both  figures  the  larger 
number  of  data  points  and  the  larger  range  of  reported  results. 

Based  upon  the  above  figures,  one  sees  that  the  reported  data  for  the  brittle  behavior  of 
out-of-plane  scaled  experiments  is  not  conclusive;  that  is,  there  is  no  clear  trend  in  the  existing 
data,  as  shown  in  Figure  2.2.  Also,  there  is  not  a  large  number  of  experiments  with  brittle 
materials,  and  no  single  experiment  showed  a  large  scale  factor.  Thus,  there  appears  to  be  room 
for  another  series  of  experiments  in  this  area. 

To  conclude  the  data  search,  we  now  turn  our  attention  to  reported  experiments  on  size 
effects  for  cyclicly  loaded  cracked  specimens;  that  is,  specimens  which  should  be  amenable  to 
reduction  by  the  Paris  Law  of  Equation  1 .4.  This  not  as  fertile  an  area  of  publication  as  in  the 
monotonic  studies.  Restricting  our  attention  to  specimen  pairs  of  identical  types  made  of  a  single 
material  with  both  width  and  thickness  scaled  in  concert,  and  which  are  tested  at  the  same 
temperature,  and  loaded  with  as  similar  frequency  and  stress  range  as  possible,  six  references 
with  admissible  data  can  be  identified  [131  to  136],  All  of  the  data  they  contain  are  for  steels  and 
the  units  used  for  C  are  those  consistent  with  expressing  a  in  inches,  and  AK  in  ksiVin.  These 


data  are  summarized  in  Table  2.1 ,  wherein  X  is  the  scaling  factor.  Clearly,  in  Table  2.1 ,  m  and  C  can 
vary  with  size:  what  is  not  clear  from  Table  2.1  are  the  implications  of  these  discrepancies  from  a 
cyclic  life  point-of-view.  To  obtain  some  idea  of  these  consequences,  we  look  to  develop  a  ratio 
for  the  lives  calculated  using  an  expression  similar  to  Equation  1 .6.  To  do  this  and  incorporate  the 
data  from  Table  2.1,  we  use  the  expression  for  lifetime,  N,  given  in  (1.5).  In  adapting  (1.5)  to 

A 

furnish  the  desired  ratio,  several  simplifying  assumptions  are  made.  First,  we  set  the  function  Y(— ) 

equal  to  one  to  facilitate  integration.  Second,  we  take  a\  as  1/32  in.,  a  common  initial  flaw  size 
taken  in  practice  because  of  detection  capabilities.  Finally  and  somewhat  arbitrarily,  we  let  a\  be 
1/4  in.,  a  value  which  can  occur  in  engineering  applications  and  one  which  allows  sufficient  crack 
growth  to  have  significant  life.  The  end  result  is; 

(  C|(m|-2)  \f mVms'2)/2-1 

N|  ~tcs(ms-2)JV  2  )  (m|-2)/2  1  ^ 

where  the  subscripts  I  and  s  refer  to  values  associated  with  the  larger  or  smaller  specimens, 
respectively,  i.e.,  ms  is  the  exponent  reported  for  the  smaller  specimen.  With  this  equation,  we 
get  a  feel  for  the  impact  of  changes  with  size. 

We  now  evaluate  (2.3)  using  the  parameters  given  in  Table  2.1  and  representative  values 
from  the  particular  sources.  We  see  that  the  resulting  ratio  may  be  greater  or  less  than  unity. 

Cases  where  Ns/N|  is  greater  than  one  may  be  interpreted  as  implying  a  conservative  estimate  of 
life  if  m,  and  C  found  via  testing  the  larger  specimen  were  used  to  predict  life  in  the  smaller,  but 
nonconservative  if  the  smaller  specimen  provided  parameters  to  predict  life  in  the  larger.  Hence 
one  immediate  consequence  is  that  neither  testing  small  and  applying  big  nor  .  ice  versa  ensures 
a  conservative  result.  To  emphasize  the  potential  for  nonconservative  estimates,  we  give  results 
from  (2.3)  as  numbers  equal  to  or  greater  than  one,  inverting  the  ratio  when  needed,  and  then 
assume  that  the  denominator  represents  the  size  used  to 
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Parameters  from  Nj/Ns;  1 .8, 1 .3,  2.4,  45.0,  1 .1 , 36.0,  1 .3 
Parameters  from  Ns/N|;  2.4,  1.1,  1.9,  1.7,  1.9,  2.6,  1.6,  2.2,  1.6 

We  can  now  grade  these  results  by  making  the  modest  requirement  that  errors  less  than  a  factor 
of  two  are  satisfactory  results,  then  regard  results  greater  than  two,  but  less  than  four,  as  poor,  and 
those  larger  than  four  as  quite  inadequate.  Based  on  this,  the  above  represents  satisfactory 
predictions  62%  of  the  time,  poor  predictions  25%,  and  inadequate  13%.  Clearly  there  is  room  for 
improvement. 

In  this  review  of  the  physical  evidence,  several  lacunas  exist  with  respect  to  our  objective 
of  gauging  how  well  the  Paris  data  reduction  scheme  performs  on  specimens  that  are  merely 
scaled.  None  of  the  cited  references  reported  using  exactly  the  same  cyclic  stress  loading  in  each 
of  their  different  sizes.  None  of  them  either  involved  sufficient  repetition  of  tests  or  presented 
results  in  such  a  way  as  to  enable  a  good  assessment  of  the  degree  of  scatter  present.  And  none 
of  them  really  allowed  the  effects  of  changing  sizes  on  the  end  points  of  the  data  fit,  AKmjn  or 
AKmax,  to  be  discerned.  Accordingly,  we  consider  a  set  of  experiments  designed  to  filling  these 
gaps,  as  well  as  those  of  the  monotonic  studies,  mentioned  earlier. 


is 


3.  MONOTONIC  LOADING  EXPERIMENTS 


In  this  chapter  we  describe  a  series  of  experiments  with  the  goal  of  showing  how  out-of- 
plane  scaling  effects  the  fracture  behavior  of  a  material.  We  begin  by  describing  a  preliminary 
experimental  study,  the  over-all  design  objectives  for  these  experiments,  and  how  they  are 
implemented.  The  discussion  then  goes  on  to  the  results  obtained  from  these  preliminary 
experiments,  and  how  these  results  lead  to  an  extended  series  of  experiments,  designed  not 
only  to  explore  out-of-plane  scaling,  but  also  to  see  if  a  correlation  might  exist  with  in-plane 
scaling.  The  chapter  concludes  with  a  presentation  of  data  obtained  from  this  expanded  series  of 
experiments. 


3.1  Preliminary  experimental  study 


The  most  important  design  objective  is  to  be  sure  that  the  experiments  do  in  fact  focus  on 
the  effects  of  out-of-plane  scaling  upon  LEFM.  Clearly  this  is  done  by  building  a  series  of 
experiments  in  which  the  only  dimension  change  is  that  of  the  thickness,  B. 


The  next  most  important  objective  is  that  the  material  used  exhibit  as  brittle  a  response  as 
possible,  since  this  response  affords  the  best  comparison  with  LEFM  by  virtue  of  complying  with 
the  underlying  assumptions  best.  However,  seeking  a  brittle  material  does  present  difficulties  in 
that  it  is  typically  costly  to  procure  in  high  quality  form,  and  hard  to  machine  into  finished 
specimens.  A  way  around  these  problems  is  to  use  a  material  that  could,  in  one  instance,  be 
ductile  for  specimen  preparation,  but  in  another  instance  be  brittle  for  testing.  One  means  to  have 
these  two  features  would  be  to  take  advantage  of  a  temperature  transition  in  a  material.  It  is  known 
that  medium  content  carbon  steels  can  behave  in  such  a  manner.  At  normal  room  temperature 
these  steels  are  ductile  and  easily  machined,  while  if  immersed  in  liquid  nitrogen  (77°  K)  they  now 
exhibit  brittle  response.  It  is  recognized  that  this  behavior  will  have  to  be  confirmed  via  a  tensile 
test  at  fracture,  and  a  description  of  this  is  below.  An  additional  advantage  to  using  this  steel  is 
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that  it  is  a  well  known  structural  material  having  a  rather  wide  range  of  applications  as  opposed  to  a 
material  like  Alumina.  The  actual  steel  chosen  is  AIS1 1045. 

Another  important  design  consideration  is  obtaining  reliable  results.  That  is  to  say.  being 
able  to  discern  real  material  response  from  experimental  noise.  Two  approaches  are  used  to 
address  this  question.  One  is  to  use  a  large  scale  factor,  changing  B  by  a  factor  of  ten  or  more.  In 
this  way  it  is  hoped  to  accentuate  any  real  physical  thickness  effects  present.  And  the  other 
approach  is  to  test  enough  specimens  so  as  to  be  able  to  estimate  in  a  statistical  fashion  the 
scatter  of  the  results  obtained. 

From  the  onset  we  also  want  these  experiments  to  follow  as  closely  as  possible  the  intent 
of  ASTM  standard  E  399  [8].  We  do  this  because  this  standard  has  achieved  quite  a  degree  of 
acceptance  with  the  fracture  mechanics  community  and  because  it  generally  represents  a  good 
method  for  obtaining  reproducible  results.  It  is  not  without  some  work  that  one  is  able  to 
incorporate  this  objective  since  the  standard  has  eleven  criteria  to  be  met.  Most  criteria  do  not 
change;  for  example:  planar  dimensions,  recording  and  reporting  requirements,  and  measuring 
methods.  None  the  less,  some  criteria  were  harder  to  keep  on  a  consistent  basis,  as  the  total 
critical  crack  length,  ac,  and  therefore  the  ratio  of  crack  length  to  specimen  width,  ac/W.  The 
criteria  that  was  impossible  to  meet  was,  of  course,  the  minimum  thickness  requirement  for  plane 
strain  testing. 

It  is  decided  that  a  disk  shaped  compact  specimen  should  be  used  for  the  fracture  tests, 
Figure  3.1 .  This  circular  specimen  allows  minimum  trim  loss  of  materials,  and  facilitates  fabrication 
procedures.  Once  a  particular  diameter  is  settled  upon,  different  thickness  are  obtained  by  simply 
slicing  off  a  length  of  stock  rod  material  roughly  equal  to  the  desired  thickness.  Having  in  mind 
testing  load  capabilities,  we  choose  a  diameter  of  2  inches.  Thereafter,  machining  is  used  to  cut 
mounting  holes,  the  initial  crack,  and  bring  the  specimens  to  final  thickness.  Eight  specimen 
thicknesses  are  produced,  from  3  inches  down  to  3/64  inches,  with  between  2  to  4  specimens  in 


each  thickness.  The  smallest  specimens  need  considerable  care  to  keep  aligned  so  that  no 
extraneous  bending  stresses  are  introduced.  Hence  a  relatively  large  sample  number,  8 
specimens,  are  prepared  to  try  to  check  that  these  stresses  were  not  present.  (A  summary 
showing  thicknesses  and  numbers  tested  is  in  Appendix  C.) 

To  gauge  the  effects  of  machining,  about  half  of  the  specimens  are  tested  in  the  "as- 
received"  state  while  the  other  half  are  heat  treated  to  reduce  residual  stress  introduced  during 
tabi,  cation.  The  heat-treating  is  such  as  to  restore  the  material  to  a  "dead  anneal"  condition  (see 
Metals  Handbook  {137],  pp  14-27).  This  is  accomplished  by  inserting  the  specimens  into  an 
argon  atmosphere  furnace  and  raising  their  temperature  to  just  above  1500°F,  holding  it  until 
equilibrium  is  reached  (one  to  two  hours  depending  upon  specimen  thickness),  and  then 
decreasing  the  temperature  at  50  degrees  per  hour  until  1200  °F,  when  the  oven  is  shut  off. 

After  heat-treatment,  all  the  specimens  are  fatigue  precracked  in  order  to  grow  a  sharp 
radiused  crack.  The  length  of  the  fatigue  precrack  is  such  that  the  ratio  of  total  crack  length,  a,  to 
width,  W,  is  within  the  desired  range,  here  0.53  to  0.55.  This  produces  adequate  crack  growth 
with  respect  to  the  ASTM  standard[8j  The  precrack  length  is  measured  using  a  traveling 
microscope.  During  this  fatigue  precracking,  the  maximum  load  is  held  within  limits  as  defined  in 
ASTM  E399  [8], 

The  actual  fracture  test  is  now  able  to  be  conducted.  The  specimens  are  mounted  in  an 
MTS  55kip  servo-hydraulic  testing  machine,  Plate  3.1 .  A  simple  liquid  nitrogen  container 
designed  to  this  end  is  then  assembled  around  the  specimen.  Once  this  is  secured,  liquid 
nitrogen  is  poured  into  the  container  until  it  covers  the  specimen,  and  is  kept  there  until  the 
specimen  reaches  a  uniform  temperature  as  based  upon  analysis  (see  Appendix  B  for 
temperature  calculations).  This  test  setup  uses  between  five  and  twenty  liters  of  liquid 
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nitrogen,  depending  upon  thickness,  to  cool  the  specimens.  The  hold  time  is  confirmed  by 
noting  the  boil  off  of  liquid  around  the  specimen.  The  testing  machine  is  then  programed  in  load 
control  to  apply  a  steadily  increasing  load  at  a  rate  such  that  the  specimen  experiences  an 
increase  of  the  stress  intensity  factor  equal  to  130  ksiVin  /min  until  fracture,  within  limits  set  in 
ASTM  E399  [8].  The  system  records  output  from  the  toad  cell  as  well  as  the  displacement,  or 
stroke,  of  the  hydraulic  ram.  The  toad  vs  stroke  curve  is  recorded. 

After  the  broken  specimen  halves  have  warmed  up,  they  are  measured  to  determine  the 
critical  crack  at  fracture.  Since  the  demarcation  between  fatigue  crack  and  precipitous  growth  to 
failure  is  usually  obvious,  it  was  fairly  easy  to  comply  with  ASTM  E399  [8]  requirements  for 
measuring  the  final  crack  front  at  the  center  and  ’quarter-thickness’  points.  The  average  of  these 
values  constitutes  ac  and  can  now  oe  combined  with  three  other  measurements,  namely 
specimen  width  (W),  thickness  (B)  and  toad  at  failure  (P*),  to  obtain  net  nominal  stress  at  fracture, 
on*.  The  net  nominal  stress  for  this  shape  specimen  is  obtained  by  remembering  that  the  stress  in 
this  remaining  ligament  results  from  both  normal  and  bending  loading  therefore  the  following 
applies:  <rn*  =  (P*/(B(W-a))  (1+  3(W+a)/(W-a)).  The  value  for  crn*  is  now  normalized  with  that  for  the 
largest  specimen  and  plotted  against  thickness  value,  also  normalized  with  respect  to  the  largest 
specimen,  it  is  interesting  to  note  that  since  the  in-plane  dimensions  are  constant,  the  variation  of 

thickness  constitutes  a  variation  of  volume  in  the  specimens.  Thus,  this  plot  of  normalized  net 
nominal  stress,  CJn*.  can  be  drawn  as  versus  thickness  or  gross  volume,  V,  where  V=  ^-BW2- 

As  stated  above,  the  brittle  nature  of  the  response  has  to  be  confirmed  by  obtaining  the 
stress-strain  curve  for  the  1045  steel  at  this  testing  temperature.  Four  specimens  are  constructed 
following  ASTM  standard  E8  {138],  see  especially  figure  7  there.  The  rectangular  cross  section  of 
these  specimens  is  0.5"  by  0.19",  and  the  distance  between  mounting  pins  is  5".  After 
manufacture,  these  tensile  specimens  are  heat-treated  in  an  identical  procedure  to  the  fracture 
specimens.  The  tensile  tests  are  carried  out  in  a  fashion  similar  to  the  liquid  nitrogen  fracture 


tests.  Results  recorded  included  the  load  displacement  curve,  yield  stress,  and  reduction  in  area 
at  fracture. 

Turning  to  the  results  obtained,  we  first  note  the  information  from  the  fracture  tests.  The 
typical  load  vs  stroke  curve  is  uniformly  linear,  after  an  initial  section  to  relieve  play  in  the  load  train, 
Figure  3.2.  This  response  is  consistent  with  desired  brittle  behavior.  Also  consistent  with  brittle 
response  is  the  absence  of  shear  lips  on  any  of  the  cooled  specimens,  all  had  essentially  flat 
fracture  surfaces.  It  is  believed  that  this  is  indicative  of  brittle  materials,  as  mentioned  in  Knott 
[139]  and  others.  Measuring  the  fractured  surfaces  showed  that  none  of  the  thickest  specimens 
have  a/W  ratios  within  the  desired  range.  In  addition,  all  specimens  are  checked  for  compliance 
with  the  crack  front  curvature  requirements  of  ASTM  standard  [8].  It  is  noted  that  crack  front  shape 
is  different  in  the  various  thicknesses,  therefore,  the  relationship  of  surface  crack  size  to  ac  varied. 

The  plot  of  normalized  net  nominal  stress  vs  normalized  gross  volume  is  shown  in  Figure 
3.3.  In  this  figure  is  shown  experimental  results  for  specimens  that  are  either  heat-treated  or  not. 
Specimens  who's  a/W  ratio  is  not  within  the  desired  range  are  identified  in  Figure  3.3.  The  curve 
of  Figure  3  3  shows  no  major  differences  between  the  behavior  of  annealed  or  non-annealed 
specimens.  Both  start  at  a  low  critical  stress,  and  as  the  volume  decreases,  the  strength 
increases.  (Details  of  results  for  all  cases  are  in  Appended  C.)  If  these  results  were  for  perfectly  or 
nearly  perfectly  brittle  response,  then  they  really  do  not  conform  with  LEFM.  This  is  because 
LEFM  attributes  thickness  effects  to  changes  in  ductility  and  otherwise  has  that  there  are  no 
effects. 


Lastly,  then,  we  consider  the  check  on  brittleness.  The  results  of  the  tensile  tests  show 
that  the  yield  strength  at  liquid  nitrogen  temperature  was  approximately  135  ksi.  However,  upon 
inspection,  the  tested  specimens  showed  excessive  amounts  of  reduction  in 
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area,  approximately  39  %.  This  would  imply  that  the  specimens  actually  lailed  in  a  ductile  fashion 
as  opposed  to  brittle.  These  results  led  to  some  questioning  of  the  material  selection  process. 
Microphotographs  of  these  tensile  specimens  show  a  marked  difference  in  the  material  from  the 
surface  towards  the  interior,  Plate  3.2.  Such  a  situation  could  come  about  as  carbon  in  the  steel  is 
lost  via  diffusion  during  the  heat-treatment.  A  loss  of  carbon  turns  a  1045  steel  into  a  1020  or 
1010  type,  having  more  ductile  properties.  As  a  result,  these  preliminary  experiments  failed  to 
meet  our  objective  of  brittle  response.  They  do,  however,  show  an  alternative  which  might  be 
successful  and  we  explore  this  next. 


3.2  Extended  experimental  study 


In  light  of  the  above  findings,  a  second  extended  series  of  experiments  are  undertaken 
with  changes  made  in  the  specimen  preparation  and  pretest  handling.  This  time  the  raw  stock  is 
annealed  first  and  then  the  specimens  fabricated  from  it.  In  this  way  the  region  containing  the 
migration  of  carbon  atoms  referred  to  above  can  be  machined  off. 


At  the  same  time  as  enacting  these  preparation  changes,  it  is  desirable  to  extend  the 
preliminary  experiments  from  simple  out-of-plane  considerations  so  as  to  explore  the  possibility  of 
a  correlation  between  in-plane  and  out-of-plane  scaling,  and  how  this  might  relate  to  the  nominal 
volume  of  the  cracked  specimen.  This  extension  of  the  experiments  can  be  carried  out  relatively 
easily  by  having  the  specimen  dimensions  change,  not  only  in  thickness  direction,  but  also  in  the 
planar,  or  width  (W)  direction.  The  selection  of  specimen  width  and  thicknesses  are  done  in  a 
fashion  so  as  to  admit  some  common  volumes  for  different  combinations  of  dimensions.  In 
particular,  for  example,  the  thickness  of  the  thinnest  large-diameter  specimen  is  chosen  such  that 
it  has  the  same  volume  as  the  second  thickest  specimen  of  the  next  smaller  diameter.  In  a 
particular  planar  size,  the  thicknesses  usually  change  by  a  factor  of  two.  In  total,  three  diameters 
were  chosen,  1  inch,  2  inches,  and  4 
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Plate  3.2  Microphotograph  (100*)  of  tensile  specimen  showing  loss  of  cartoon  rich  phase 
near  surface. 


inches;  thicknesses  ranged  from  1  and  1/2  inches  down  to  3/64  inch.  The  actual  sizes  of 
specimens  fabricated  are  shown  as  a  matrix  in  Figure  3.4,  eight  gross  volumes  are  tested. 

To  start  this  extended  experimental  activity,  a  bulk  lot  of  1045  steel  was  purchased.  It 
came  as  a  4-inch  diameter  rod.  This  rod  was  cut  up  into  1-foot  lengths  for  heat-treating.  After 
these  1-foot  sections  are  heat-  treated  to  a  "dead  anneal"  condition,  as  in  the  preliminary 
experiments,  the  4-inch  diameter  is  then  reduced  to  either  1-inch  or  2-inch  as  needed,  or  if  a  4- 
inch  specimen  is  required,  the  diameter  is  simply  trued.  Although  time-consuming,  this  method  of 
fabrication  has  the  advantage  that  all  specimens  were  made  from  the  same  stock  of  material,  and 
also  that  the  crack  front  in  all  specimens  was  at  about  the  same  location  with  respect  to  the  original 
4-inch  diameter  stock.  After  the  proper  specimen  diameters  were  obtained,  the  specimen 
preparation  was  handled  in  a  similar  fashion  to  that  of  the  preliminary  experiments  using  scaied 
slitting  saws  as  needed.  The  final  fabrication  step  is  grinding  to  final  thickness. 

Now,  the  specimens  are  fatigue-cracked  to  the  same  a/W  value  as  the  preliminary  series. 
Fatigue  pre-cracking  proceeds  as  in  the  preliminary  experiments  except  that  it  becomes  more 
difficult  to  meet  all  the  criteria  of  ASTM  E399  since  some  are  expressed  in  terms  of  absolute 
values;  to  wit,  the  minimum  precrack  is  0.05  inches  [8],  A  more  important  problem  with  the 
standard  [8],  is  the  application  of  AK  to  the  precrack  procedure  involving  various  specimen  sizes. 
When  two  specimens  with  different  W's  are  precracked  with  the  same  AK,  the  smaller  specimen 
experiences  a  larger  AOn,  since  AOn  OC  AK/Vw  for  constant  a/W.  In  an  effort  to  minimize  the 

strain-hardening  experienced  by  the  specimens  before  the  actual  fracture  test,  it  is  decided  to 
use  a  minimum  AOn  in  the  crack  growth  procedure.  All  precracking  was  started  at  On  =  27  ksi  and 

slowly  increased  until  97%  of  the  desired  crack  growth  could  occur  within  several  hundred 
thousand  cycles.  For  the  last  2  to  3%  of  crack  growth  the  On  was  reduced,  meeting  the  intent  of 

the  ASTM  standard  [8j. 
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The  order  of  fracture  testing  of  the  extended  series  of  specimens  is  done  in  a  somewhat 
random  fashion  so  that  trends  in  the  data  will  not  be  associated  with  testing  order.  First  a  single 
specimen  from  a  cross  section  of  sizes  is  tested.  Then,  the  testing  can  be  thought  of  as  going 
along  by  gross  volumes.  The  eight  gross  volumes  are  collected  into  3  groups  by  size,  labeled 
large,  middle  and  small,  with  3,  3  and  2,  respectively,  volumes  in  the  groups.  The  order  of  group 
testing  then  is:  middle,  large,  small,  middle,  large,  and  finally,  small.  The  actual  liquid  nitrogen 
fracture  testing  is  done  in  a  similar  manner  to  that  of  the  preliminary  tests,  as  are  the  post-fracture 
measurements. 


In  addition,  with  this  new  stock,  we  repeat  the  tensile  tests  at  liquid  nitrogen 
temperatures.  To  accomplish  this  a  heat-treated  section  is  quartered  and  tensile  specimens 
prepared  along  the  lines  of  ASTM  A  370  [140],  Four  specimens,  having  diameters  of  either  0.35 
or  0.175  inches,  are  prepared.  They  are  tested  in  a  liquid  nitrogen  bath,  just  as  those  of  the 
preliminary  series. 


We  are  now  in  a  position  to  present  the  results  obtained  for  the  extended  experimental 
study,  starting  with  the  tensile  tests.  This  time  the  specimens  have  between  4.6%  and  2.6% 
reduction  in  area  and  ultimate  strengths  of  158.0  to  161 .6  ksi.  A  typical  stress-strain  cun/e  is 
shown  in  Figure  3.5.  These  results  are  comparable  with  data  from  MacGregor  [141]  for  the  same 
material,  and  seem  to  represent  very  brittle  response. 


The  fractured  specimens  of  this  series,  in  general,  are  redolent  of  those  of  the  preliminary 
study,  having  no  shear  lips  and  relatively  flat  fracture  planes,  with  patent  demarcation  between 
areas  of  fatigue  or  precipitous  cracking,  see  Plate  3.3.  A  typical  load-displacement  curve  for  a 
specimen  in  this  series  is  rather  similar  to  one  of  the  preliminary  study. 
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Upon  review  of  the  results  from  this  series,  it  is  discovered  that  the  curvature  of  the  critical 
crack  is  more  than  anticipated,  although  not  always  enough  to  invalidate  the  test  with  respect  to 
the  standard  [8].  As  a  result,  many  of  the  final  raw  a/W  values  are  not  within  the  closely  confined 
region  desired.  To  overcome  this  problem,  it  was  decided  that  a  weighted  average  be  used  with 
the  results. 

Upon  examining  the  data  further,  a  mean  value  of  a/W  of  0.553  seemed  to  include  as 
many  tests  as  possible.  The  actual  procedure  for  weighting  test  results  to  arrive  at  this  mean  for 
each  size  is  as  follows.  First  the  participation  factor,  f,is  determined.  This  is  done  by  setting  f  =1 
when  0.540  <  a/W  <  0.566,  0  <  f  <  1  when  0.466  <  a/W  <  0.540  or  0.566  <  a/W  <  0.640  and  such 
that  the  weighted  mean  equals  0.553,  and  f  =  0  for  other  a/W.  Under  these  rules  there  still 
remains  some  flexibility  in  choosing  f ;  to  reduce  this  we  tried  to  pick  f  such  that  the  maximum 
participation  of  test  data  occurred.  The  actual  choice  of  fs  lead  to  a  sum  of  f’s  equal  to  30.46,  with 
some  participation  from  42  out  of  53  specimens  tested,  and  all  sizes  represented.  With  these  fs, 
the  average  response  for  the  specimen  size  is  computed  as  the  sum  of  the  products  of  on*  and  f 
for  each  specimen,  divided  by  the  sum  of  the  fs.  (Details  of  results  and  participation  factors  are  in 
Appendix  C.) 

For  clarity,  we  first  present  the  results  in  terms  of  thickness,  or  out-of-plane  scaling. 

Figure  3.6  is  a  plot  of  normalized  strength,  on*,  as  in  Figure  3.3,  for  each  specimen  of  a  given 

diameter,  versus  normalized  thickness.  The  normalizing  factors  come  from  values  of  the  largest 
thickness  in  each  diameter.  If  the  information  of  this  figure  were  to  follow  the  LEFM  thickness 
prediction  discussed  in  chapter  1,  all  the  points  should  be  along  a  horizontal  line  through  1 .0,  but 
they  are  not. 

We  now  consider  in-plane  and  out-of-plane  scaling  at  once.  To  enable  in-plane  and  out- 
of-plane  size  effects  on  strength,  Gn*,  to  be  plotted  on  a  single  graph,  we  use  gross  volume,  as 

before.  That  is,  we  plot  On*  versus  V  as  seen  in  Figure  3.7.  As  with  the  preliminary  study,  the 
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Plot  of  normalized  net  strength  vs  normalized  thickness, 
extended  study. 


normalizing  factor  come  from  the  largest  specimen.  Included  with  the  stress  values  for  eacn 
volume,  are  range  bars  for  the  results.  In  general,  this  figure  portrays  the  type  of  behavior 
expected:  that  is  as  the  specimen  size  gets  small,  the  strength  increases.  This  decrease  in  size  is 
due  to  both  in-plane  and  out-of -plane  changes,  the  in-plane  changes  having  the  more 
pronounced  effect. 

Also  shown  in  Figure  3.7,  as  the  dotted  line,  is  the  prediction  of  strength  based  upon 
LEFM.  The  prediction  uses  the  largest  specimen  as  the  starting  point.  The  lines  are  flat  since 
there  is  no  thickness  effected  predicted.  Note  that  as  the  specimens  get  smaller,  the  predictions 
get  further  off  the  mark,  missing  five  of  the  eleven  ranges  entirely,  and  only  getting  the  ranges 
for  the  smallest  specimens  because  of  the  large  amount  of  scatter  involved.  The  mean  error  for 
the  eleven  predictions  is  7%  while  the  largest  single  error  is  an  over  prediction  of  about  16%.  We 
can  get  a  feeling  for  how  well  LEFM  does  with  predicting  the  thickness  effect  alone  by  assuming 
LEFM  predicts  the  On*  of  the  largest  gross  volume  in  each  planar  set  and  looking  at  the  resulting 

errors.  The  mean  error  for  eight  predictions  is  10%  and  the  largest  single  error  is  an  under 
prediction  of  19%. 

The  dashed  line  of  Figure  3.7  represents  a  Weibull  type  plot  [142],  where  the  normalized 
stress  is  a  function  of  the  volume  raised  to  some  value,  i.e.  On*=  (V)a.  The  value  of  the 

exponent,  a,  is  selected  as  the  one  that  best  fits  the  experimental  data  in  a  least-squares  fashion, 
weighted  by  the  participation  factor,  here  a  =  -0.1 183  .  This  method  does  about  the  same  as 
LEFM  in  that  it  misses  the  same  number  of  ranges,  and  the  average  error  is  0.08,  and  the  worst 
single  error  is  an  under  prediction  of  20%. 

Finally,  the  data  is  fitted  to  another  model,  as  shown  in  the  solid  curve  of  Figure  3.7.  This 

model  represents  the  normalized  strength  as  a  function  of  a  product  of  the  in-plane  scaling  factor 
raised  to  a  power  and  the  out-of-plane  factor  raised  to  another;  i.e.  On’  =  (W/Wo)“  (B/Bo)P, 

where  the  0  subscript  indicates  the  normalizing  value,  or  the  value  of  the  largest  gross  volume. 


Plot  of  normalized  net  strength  vs  normalized  gross  volume, 
extended  study. 


The  two  exponents,  a  and  (3,  are  selected  so  as  to  hit  as  many  data  ranges  as  possible,  here  ten 
out  of  the  eleven.  Of  the  two  exponents,  a  is  largest,  -0.31  versus  -0.04  for  p,  and  therefore  has 
the  larger  impact  on  the  strength.  With  these  values,  this  model  has  a  mean  error  of  5%,  with  the 
worst  single  error  being  a  1 6%  under  prediction.  It  is  interesting  to  note  that  these  two  values  are 
in  the  same  trend  as  LEFM:  that  is  to  say  -0.3  is  not  far  from  -0.5,  and  -0.04  is  close  to  zero.  A 
possible  explanation  for  the  success  of  this  model  is  that  the  surfaces  may  be  different  from  the 
interior  by  virtue  of  the  machining  they  obtained. 


4.  CYCLIC  LOADING  EXPERIMENTS 

Herein  we  describe  a  series  of  experiments  whose  primary  objective  is  to  be  a  pilot  study 
on  scaling  effects  in  cyclic  loading  situations.  We  start  this  chapter  by  describing  the  experimental 
objectives.  This  is  followed  by  a  description  ot  how  the  specimens  are  prepared  and  the 
experiments  conducted.  We  conclude  the  chapter  by  discussing  the  data  reduction  methods 
used  and  results  found. 

The  primary  objective  of  this  experimental  series  is  to  conduct  a  pilot  study  into  the  scaling 
effects  in  a  cyclic  loading  situation,  and  see  how  consistent  these  results  are  with  LEFM.  We 
chose  to  look  at  the  complete  scaled  situation,  i.e.  proportional  changes  in  length  and  thickness 
as  well  as  width,  as  shown  in  Figure  1 .3.  We  attempt  to  isolate  other  factors  by  doing  the 
following:  using  the  same  material,  cycling  at  the  same  frequency,  conducting  the  tests  at  the 
same  temperature,  applying  a  constant  stress  amplitude  in  all  tests,  and  employing  the  same  R 
ratio,  where  R  is  the  minimum  cyclic  stress  divided  by  maximum  cyclic  stress.  As  with  the 
monotonic  experiments,  we  wish  to  control  scatter  and  therefore  will  use  several  specimens  in 
each  size.  Also,  we  attempt  to  comply  with  ASTM  standards,  namely  E  647  [143]  because  of  their 
general  acceptance  within  the  fracture  mechanics  community. 

The  experiments  use  center-cracked  panels  scaled  with  respect  to  width,  thickness  and 
length  by  a  factor  of  four,  a  factor  of  four  being  the  largest  compatible  with  the  test  rig  and 
machining  practices,  Figure  4.1.  Detailed  preparation  of  the  specimens  is  in  accordance  with 
ASTM  Standard  E  647  [143].  The  center  starter  crack  is  machined  into  the  material  using  a 
wedge-shaped  flycutter.  The  result  of  this  is  that  the  both  ends  of  the  crack  have  a  chevron 
starter  area.  This  method  of  producing  the  center  crack  is  scaled  for  both  sizes,  i.e.,  an  exactly 
1  /4-size  flycutter  is  used  to  cut  the  smaller  specimen.  During  the  fabrication,  the  specimens  are  all 
aligned  in  such  a  manner  that  the  rolling  direction  is  the  same  with  respect  to  the  center  crack. 

The  next  step  in  the  preparation  is  to  measure  the  specimens  and  the  center-crack  geometry 
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completely.  Inspect’on  microscopes  are  used  for  this  operation.  The  specimens  are  made  ot 
aluminum  2024-T3  in  contrast  to  the  steel  of  the  monotonic  loading  experiments.  Loading  for 

both  specimen  sizes  is  at  the  same  frequency  (20  Hz),  with  an  R-ratio  of  0.1,  and  applied  stress 
variation,  AG,  equal  to  9.6  ksi.  Testing  is  carried  out  on  sets  of  six  specimens  of  each  size,  six 

being  the  minimum  deemed  necessary  to  gauge  scatter. 

The  execution  of  the  experiments  tried  to  follow  the  guidelines  set  forth  in  ASTM  E  647 
[143].  The  experiments  are  conducted  on  a  typical  MTS  22  Kip  servo-hydraulic  feedback  testing 
machine.  The  test  set-up  is  shown  in  Plate  4.1.  The  order  of  testing  is  such  that  some  large 
specimens  are  tested,  then  smaller  specimens,  finally  finishing  with  the  remainder  of  the  large 
specimens.  This  is  done  to  limit  any  systemic  trend  in  the  data  due  to  testing  sequence.  Another 
aspect  of  the  experimental  procedure  is  the  measurement  of  the  crack  during  cycling.  To 
accomplish  this,  an  optical  method  is  used.  This  consists  of  a  strobe  light  synchronized  to  the 
maximum  load  signal  of  the  machine  and  a  traveling  microscope.  The  operational  details  of  this 
setup  are  that  the  measurements  are  taken  only  from  the  front  surface  and  that  the  left  and  right 
sides  of  the  crack  are  alternately  measured.  Since  a  strobe  is  used,  there  is  no  need  to  halt  the 
machine  operation  for  measurement.  The  experiments  then  continued  until  fracture.  Cycles  to 
fracture,  N,  are  recorded  as  well  as  the  crack  length  before  catastrophic  failure,  af .  (A  summary  of 
the  raw  test  data  for  all  specimens  is  contained  in  Appendix  C.) 

The  next  step  is  to  do  the  data  reduction.  This  starts  by  combining  the  left  and  right  side 
measurements  for  a  specimen.  Incorporated  in  this  combining  is  the  implicit  assumption  that  crack 
growth  is  symmetric,  being  hopefully  justified  by  the  efforts  to  ensure  symmetry.  The  ASTM 
standard  [143]  presents  two  possible  ways  for  data  reduction.  These  two  methods  are  the  secant 
method  and  the  seven  -point  polynomial  method.  Both  were  applied  to  the  data. 

Individual  results  in  the  form  of  points  on  a  da/dN  versus  AK  plot  are  shown  in  Figure  4.2. 
This  figure  superimposes  the  results  for  all  specimens  in  the  two  sizes.  Viewing  this  figure,  one 


notes  the  typical  three-stage  shape  of  the  two  superimposed  curves.  Although  the  center  linear 
portions  of  the  two  curves  are  roughly  parallel,  it  is  important  to  note  that  »he  initial  and  final  stages 
of  the  two  curves  have  a  relatively  large  offset.  The  Paris  parameters  of  {1 .4),  obtained  fr  m  a  lest- 
squares  fit  to  the  data,  are  shown  in  Table  4.1.  Both  Figure  4.2  and  the  results  in  Table  4.1 
essentially  remain  the  same  whether  the  secant  method  or  the  seven-point  polynomial  method  is 
used. 


The  first  thing  to  note  concerning  the  values  in  Table  4.1 ,  is  that  it  would  not  appear 
reasonable  to  attribute  the  changes  in  m,  or  C  with  size  to  scatter  alone,  since  neither  the  mean  m 
nor  the  mean  C  for  one  size  lies  within  the  corresponding  range  for  m,  or  C  of  the  other.  Indeed, 
the  ranges  themselves  for  m  and  C  overlap  less  than  16%  and  24%  of  their  combined  extents, 
respectively.  It  thus  seems  likely  that  these  results  demonstrate  that  m  and  C  can  exhibit 
significant  dependence  on  size.  Furthermore,  the  end  points  of  the  validity  of  the  data  fits  are 
size  dependent.  Hence  if  one  were  estimating  what  would  occur  in  a  small  specimen  using  the 
data  from  the  big  when  crack  propagation  started  at  a  AK  of,  say,  13,  one  one  might  be  inclined  to 
predict  indefinite  life,  or  at  least  to  predict  it  lasting  as  many  cycles  as  the  larger  specimens  typically 
survived  -  here,  on  average,  27,000.  In  fact,  the  smaller  specimens  only  lasted  an  average  of 
4,000  cycles  from  this  point.  While  it  is  now  generally  recognized  that  the  notion  of  a  threshold  AK 
as  a  material  property  is  invalid,  this  example  does  serve  notice  of  the  dangers  of  regarding  either 
AKmin  or  AKmax  as  a  size-independent  material  constant  .  In  all,  then,  for  the  present  set  of 
experiments,  cyclic  life  predictions  based  on  the  Paris  data  reduction  scheme  may  be 
unacceptably  unreliable. 


t  Such  a  potentially  dangerous  interpretation  is  implicit  in  the  sinh  fit  of  Annis,  Wallace, 
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and  Sims  [144], 
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Figure  4.2  Crack  growth  rates  versus  AK. 


Finally,  here  we  check  the  size  dependence  given  earlier,  namely, 


N_s  =  x(nV2  -1) 


We  evaluate  (4.1)  with  the  data  of  Table  4.1 ,  using  m  from  both  sizes.  The  predicted  value  of 
(4.1),  using  m  from  the  small  and  large  specimens  is  1 .66  and  2.36,  respectively.  The  ratio  of 
actual  average  N*s  to  actual  average  N*|  is  2.93.  So  neither  represents  very  good  agreement. 


Quantity 


m,  mean 

m,  range 

C,  mean 
(xIO10) 

C, range 
(xIO10) 

N,  mean 
(xIO3) 

AKmjn 

(ksiVin) 

AKmax 

(ksiVin) 


Small 

specimens 

(fc.1) 


2.59-3.03 


91-228 


Large 

Specimens 

(k=4) 


2.88-3.54 


20-140 


5.  CONCLUDING  REMARKS 


In  this  final  chapter,  we  attempt  to  summarize  the  information  contained  in  the  preceding 
chapters.and  discuss  its  implications  for  LEFM.  We  begin  by  looking  at  the  three  particular  scaling 
situations  in  question  and  review  the  information  obtained  from  both  literature  and  experimental 
studies  in  each  situation.  We  continue  by  expressing  some  thoughts  on  the  consequences  of 
this  information  for  LEFM,  and  conclude  by  discussing  possible  areas  of  further  activities. 

We  first  address  the  situation  of  in-plane  scaling  for  monotonically-loaded  specimens.  As 
expressed  in  Chapter  2,  the  summary  of  literature  review  by  Sinclair  and  Chambers  [7]  indicates 
that  reported  data  do  not  agree  well  with  the  square-root  scaling  prediction  of  LEFM  see  (1.3),  as 
shown  in  Figure  2.1 .  The  prediction  becomes  further  removed  from  the  data  as  the  scaling  factor 
increases,  tending  to  over  predict  the  strength  in  smaller  specimens.  The  experimental  results 
described  in  Chapter  3  also  showed  a  failure  to  predict  strength  changes  with  scaling.  However, 
in  these  experiments,  the  overestimate  occurred  as  the  specimens  got  smaller,  and  the 
inaccuracies  were  not  as  pronounced  as  for  some  of  the  data  found  in  [7].  Nonetheless,  in  all 
LEFM's  performance  in  this  regard  is  less  than  satisfactory. 

Next,  we  turn  our  attention  to  the  out-of-plane  scaled  situation.  The  literature  review  in 
this  area  indicates  that  there  is  no  definitive  trend  to  the  reported  data.  In  fairness,  it  should  be 
pointed  out  that  these  data  represent  a  limited  number  of  experiments,  at  least  for  brittle  materials, 
and  which  typically  did  not  have  large  scaling  factors.  The  experiments  reported  upon  in  Chapter 
3  exhibited  a  strength  variation  with  out-of-plane  scaling.  Figure  3.7  shows  that  the  out-of-plane 
scaling  effect  does  not  have  as  pronounced  an  impact  on  the  normalized  strength  as  does  the  in¬ 
plane  effects.  For  this  case,  the  "flat  line"  prediction  of  LEFM  is  conservative  and  underestimates 
the  strength  increase.  However,  because  of  the  limited  ductility  of  the  material  in  the  test 
conditions,  it  is  clear  that  no  explanation  which  is  based  solely  on  a  material  ductility  aspect  could 
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fully  explain  this  situation.  In  all,  the  understanding  of  thickness  effects  in  LEFM  is  less  than 
complete. 


To  quantify  the  above  situation  for  in-plane  and  out-of-plane  scaling,  consider  the 
expression  on*  =  Xa,  and  the  changes  in  On*  as  a  result  of  setting  a  to  either  0.3  or  0.5.  We 
construct  a  ratio  of  A.0-3/  A.0-5,  and  note  how  this  ratio  varies  with  the  scaling  factor,  X.  For  X  =  2,  the 
ratio  is  1 .23/1 .41 ,  or  0.87,  indicating  about  a  13%  difference  in  predicted  strength.  As  X  is 
increased,  the  ratio  changes,  such  that  when  X  =  8,  the  ratio  is  1 .87/2.83,  or  0.66.  If  X  gets  up  to 
16,  the  ratio  becomes  0.57.  Clearly  the  impact  of  a  is  dependent  upon  the  scaling  factor;  if  X  is 
large,  large  differences  in  strength  predictions  result  from  small  variations  in  X.  However,  for  out- 
of-plane  scaling,  with  a  being  eitherO.O  or  0.04,  the  above  ratio  will  only  vary  by  12%  even  for  a 
scaling  of  20. 

Finally  we  address  the  last  question  under  study,  we  consider  scaling  effects  in  cyclically- 
loaded  situations.  The  literature  review  in  this  area  was  most  noteworthy  for  its  relatively  few 
reports.  These  contributions  seem  to  question  the  validity  of  assuming  the  invariability  of  the 
material  parameters,  C  and  m,  of  the  Paris  Law,  (1.4).  In  using  such  reports  to  determine  the  effect 
of  scaling  on  lifetime  predictions,  there  seems  to  be  no  clear  trend.  Sometimes  the  data  indicates 
increases  in  lifetime,  while,  in  other  cases,  decreases  in  lifetime.  Also  noteworthy  in  these 
contributions  is  the  lack  of  information  about  behavior  of  the  range  of  applicability  with  scaling. 

The  reported  experimental  results  of  Chapter  4  raised  similar  questions  about  C  and  m  and  the 
ability  to  predict  lifetimes.  However,  what  may  be  of  more  important  consideration  is  that  these 
experiments  indicated  a  possible  shift  of  the  range  of  applicability.  This  could  imply  a  finite  life  for 
items  believed  to  be  designed  within  "infinite  life"  parameters.  In  sum,  the  LEFM  treatment  of 
cyclic  life  prediction  for  scaled  specimens  appears  to  be  potentially  quite  unreliable. 

Overall,  LEFM  is  trendwise  generally  correct  in  its  predictions  of  size  effects,  but  so 
oversimplified  that  it  is  unreliable  as  an  engineering  technology.  It  seems  that  what  ever  size 
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effects  are  present  for  a  given  material,  LEFM  does  not  usually  predict  them  accurately  enough. 
This  is  probably  because  these  effects  are  dependent  upon  the  microstructure  of  the  material 
under  consideration  and  LEFM,  being  based  on  the  continuum  mechanics  of  homogeneous 
materials,  does  not  include  the  effects  of  microstructure.  In  order  to  accommodate  size  effects, 
therefore,  microstructure  will  need  to  be  incorporated.  Just  how  this  should  be  done  is  a 
challenging  question,  one  which  is  beyond  the  scope  of  the  present  work.  However,  it  seems 
fairly  certain  that  such  a  question  will  have  to  be  considered  by  future  workers  if  we  are  to  arrive  at  a 
more  reliable  technology. 
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APPENDIX  A.  TABLES  OF  THICKNESS  DATA 

The  thickness  data  gleaned  from  the  references  is  herein  presented  in  a  format  to 
complement  other  local  work  [1.7].  The  following  tables  present  data  grouped  with  respect  to 
material  tested,  i.e.,  steels,  aluminum  alloys,  other  metals,  and  non-metals.  In  each  of  these, 
the  table  format  is  the  same  and  begins  with  the  contributor's  surname  and  reference  number 
Below  this  is  a  generic  identification  of  the  specific  material  reported,  e  g.,  for  steel,  4340; 
aluminum,  2024-T635;  for  other  metals,  Tj6A-6V-2Sn,  etc..  These  entries  are  followed  by  an 

indication  of  specimen  used  according  to  the  following  abbreviations:  CCP.. .center  cracked 
panel,  CTS... compact  tension  specimen,  RCT... found  compact  tension,  RDS... rotating  disk 
specimen,  SEN... single-edge  notch,  WOL.. .wedge  opening  loading,  3PB... three-point  bend, 

4PB... four-point  bend. 

The  number  in  parentheses  following  one  of  the  above  abbreviations  indicates  the 
total  number  of  specimens  reported.  If  a  further  number  follows,  it  represents  the  testing 
temperature:  no  further  such  number  indicates  tests  performed  at  room  temperature.  The 
fourth  column  in  the  table  represents  the  normalizing  or  representative  thickness,  B0.  The 
next  column  is  the  thickness  scaling  factor  for  the  specimens  reported,  Xg.  The  fourth 
column  is  the  ratio  of  reported  critical  strength  to  that  at  a  thickness  of  B0,  i.e.,  Oj*/Oo* 

The  final  column  classifies  the  response  in  accordance  with:  b... brittle,  b-d... brittle  ductile,  and 

::1 


d... ductile. 


A, 


Source,  material,  specimen 
type  (#  tested) 


Andrews  et  al.  [9], 

A469,CD[2-6] 

1.0 

0.38 

0.77 

1.54 

1.35 

0.55 

1.11 

2.22 

Bandyopadhyay  et  al.  [1 OJ, 

0.22  C,  0.87  Mn,  0.07  Si,  0.02  S,  &  0.02P; 

3PB  [18] 

0.045 

5.00 

2.73 

1.32 

8.70 

22.12 

0.071 

1.67 

3.17 

5.71 

0.151 

0.75 

1.46 

6.49 

0.215 

1.11 

1  82 

0.245 

0.47 

0.96 

1.55 

4.32 

Banerjee  [11], 

ASTM  516-Gr70,  CT 

0.5 

0.20 

0.39 

0.67 

1.0 

2.0 

Barnby  and  Al-Daimalani  [12], 

BS1956A  C-ICr,  3PB  (10  -  30) 

0.44 

0.89 

1.33 

1.78 

2.22 

2.67 

3.56 

BS1348E 

1/2Cr-1/2Mo-1/2V,3Pb  (10  -30) 

0.44 

0.89 

1.33 

lication 


BS1458A 

I^Mn-Ni-Gr-Mo,  3PB  (10  -30) 


BS1458B 

11/2Ni-Cr-Mo,  3PB(15&25) 

Barnby  and  Al-Daimalani  [13], 
BS1760B,  3PB  (25  &  30) 


BS1456A 

11/2Mn,3PB(10&30) 

Batte  et  al  [15], 

1%Cr  Mo  VIP,  SEN 

Brown  and  Sprawley  [18], 
Maraging  steel, 4PB 


Maraging  steel, SEC 

Maraging  steel, CCP 

Maraging  steel, 4PB 
Maraging  steel,  SEC 

Maraging  steel,  CCP 

Maraging  steel,  4PB 


Table  A.l  (con't) 


Gj 

Source,  material,  specimen 

Temp. 

Bo 

*B 

Class 

GO* 

3.64 

0.97 

b 

Chell  and  Worthington  [20], 

1.37 

0.36 

0.99 

d 

Ducol  W30A,  CTS,  (8) 

0.73 

0.95 

d 

1.09 

1.01 

d 

1.45 

0.90 

d 

1.00 

0.50 

1.13 

d 

1.00 

1.00 

d 

1.50 

1.06 

d 

2.00 

0.92 

d 

Elsender  et  al.  [21], 

0.45 

8.69 

1.20 

d 

1%CrMoVlP,  SENT 

16.95 

1.27 

d 

Fukakura  [23], 

-148°F 

0.10 

0.80 

1.09 

d 

0.15  C,  0.01  Si,  0.71  Mn,  0.01  P, 

2.00 

0.77 

d 

0.02  S;  4PB 

4.00 

0.70 

b-d 

6.00 

0.71 

b-d 

8.00 

0.78 

d 

12.00 

0.72 

d 

Griffis  [25],* 

0.125 

1.91 

1.08 

d 

HY-80,  3PB 

0.97 

1.0 

d 

0.49 

0.86 

d 

Hollstein  et  al.  [27], 

0.29 

2.72 

1.02 

d 

St  E  460,  SEN 

5.44 

1.04 

d 

0.47 

0.86 

0.99 

d 

1.72 

1.03 

d 

3.43 

1.08 

d 

0  63 

0  62 

0.98 

d 

1.23 

1.01 

d 

2.46 

1.01 

d 

0.79 

0.50 

1.03 

d 

0.99 

1.00 

d 

1.98 

0.97 

d 

0.94 

0.42 

0.96 

d 

0.84 

0.99 

d 

Huang  and  Gelles  [28], 

0.25 

0.39 

1.00 

d 

HT-9,  CT 

1.19 

1.00 

d 

1.87 

1.03 

d 

Jones  and  Brown  [32] 


0.135  0.37 


0.94 


b-d 


Table  A.1  (con’t) 


Oj 

Source,  material,  specimen 

Temp. 

Bo 

Class 

O0‘ 

-184°F 

1.2 

1.0 

d 

2.4 

1.02 

d 

3.6 

0.60 

d 

4.8 

0.57 

d 

6.0 

0.96 

d 

Mowbray  et  al.  [42], 

~-50°F 

0.04 

5.0 

1.23 

d 

A-3023,  3PC 

10.0 

1.40 

d 

-100°F 

0.04 

5.0 

1.05 

d 

10.0 

1.11 

d 

-150°F 

0.04 

5.0 

0.84 

d 

10.0 

0.73 

d 

Nl-Mo-V,  SEN 

-100°F 

0.164 

0.12 

1.75 

d 

0.76 

1.09 

d 

-320° F 

0.164 

0.12 

1.7 

d 

0.76 

1.22 

d 

Munz  and  Keller  [43], 

-184°F 

0.6 

0.20 

0.93 

b-d 

35Ni-CR-Mo-16 

0.40 

1.02 

b-d 

0.67 

1.02 

b-d 

1.00 

1.00 

b-d 

1.67 

1.02 

b-d 

-166°F 

0.6 

0.20 

0.88 

b-d 

1.00 

1.00 

b-d 

1.67 

1.00 

b-d 

-148°F 

0.6 

0.20 

0.72 

b-d 

0.40 

0.76 

b-d 

0.67 

0.96 

b-d 

1.00 

1.00 

b-d 

-130°F 

0.6 

0.20 

0.72 

b-d 

0.40 

1.05 

b-d 

0.66 

0.95 

b-d 

1.00 

1.00 

b-d 

1.66 

1.28 

b-d 

-112°F 

0.6 

0.66 

1.06 

b-d 

1.00 

1.00 

b-d 

1.66 

1.25 

b-d 

Neale  [44], 

1.03 

1.89 

0.98 

b-d 

Table  A.i  (con't) 


Source,  material,  specimen 

Temp. 

Bo 

0.50 

1.00 

2.00 

300M 

0.125 

2.00 

4.00 

0.25 

1.00 

2.00 

4.00 

0.50 

1.00 

2.00 

Sliney,  Jr.  [55], 

5%  Cr-Mo-V,  DEN 

200°F 

0.075 

0.53 

0.80 

1.07 

1.33 

76°F 

0.075 

0.27 

0.53 

0.80 

1.07 

1.33 

-40°F 

0.075 

0.27 

0.53 

1.07 

1.33 

-65°F 

0.075 

0.27 

0.53 

0.08 

-100°F 

0.075 

0.27 

0.53 

CCP 

76°  F 

0.857 

0.05 

0.14 

0.19 

0.25 

Srawley  and  Beachem  [57], 

G8  steel,  CCP 

80°F 

0.25 

0.52 

0.40 

0.32 

0.24 

0.16 

0.08 
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Table  A.1  (con't) 


Source,  material,  specimen 


68  Steel,  CCP 


W2,  CCP 


Steigerwald  [58], 
H-11,  CCP 


Steigerwald  and  Hanna  [59], 
AM  355  SST,  CCP 


Sullivan  and  Stoop  [60], 
4130,  CCP 


Temp. 

Bo 

0°F 

0.25 

0.40 

0.20 

-100°F 

0.25 

0.40 

0.20 

-200°F 

0.25 

0.40 

0.20 

-280°F 

0.25 

0.40 

0.20 

+300° F 

0.25 

0.40 

0.32 

0.24 

0.16 

0.08 

80°  F 

0.25 

0.40 

0.32 

0.24 

0.16 

0.08 

-280° F 

0.25 

0.40 

0.32 

0.24 

0.16 

0.08 

0.175 

0.09 

0.15 

0.21 

0.30 

0.55 

0.175 

0.01 

0.01 

0.02 

0.02 

0.03 

0.97 

0.03 

0.06 

1.68 

Table  A.1  (con't) 


Source,  material,  specimen 


D6A,  CCP 


Taira  and  Tanaka  [62], 
Cr-Mo-V,  CT 


Winne  and  Wundt  [65], 
Ni-Mo-V,  RDS 

Ni-Mo-V  (0.27  Mo),  RDS 
Cr-Mo-V 


Yukawa  [66], 
Alloy  steel 


Cl  Cl  Q.  Cl  CL  O.  Cl  Cl  Q.  Q.  Cl  Q. 


Table  A.2  (con't) 


Source,  material,  specimen 


2024-T351,  CCP 


Frediani  [22], 
2219-T851,  CCP 


Grewal  and  Weiss  [24], 
2024-T351 ,  SEN 


Heyer  and  McCabe  [26], 
7075-T6,  CCP 

Hudson  and  Lewis  [29], 
2219-T851 ,  Not  stated 


Bo 

gi 

oo‘ 

0.65 

0.06 

1.06 

0.12 

1.17 

0.18 

1.18 

0.24 

1.18 

0.30 

1.24 

0.36 

1.16 

0.48 

1.09 

0.61 

1.06 

0.22 

0.45 

1.01 

0.72 

1.00 

0.36 

0.13 

0.88 

0.27 

0.89 

0.43 

0.93 

0.43 

0.11 

1.04 

0  23 

1.03 

0.57 

0.08 

0.91 

0.17 

0.89 

0.28 

0.92 

0  64 

0.07 

0.95 

0.15 

1.00 

0.25 

1.00 

0.78 

0.06 

0.78 

0.13 

0.78 

0.20 

0.82 

0.25 

1.00 

1.00 

0.50 

1.07 

0.25 

1.15 

0.25 

1.00 

1.00 

0.50 

0.91 

0.25 

0.82 

0.38 

0.17 

1.55 

0.08 

1.84 

1.50 

1.00 

1.02 

0.83 

1.63 

0.58 

1.98 

a  a.  <x  o.  a.  a.  a.  a  a  a  a  cl  a.  a  o.  a.  a.  a.  a  a.  cl  cl  q.  a.  a.  o.  cl  aaaaaaaa 


Table  A.2  (con't) 

Source,  material,  specimen 

Irwin  et  al.  [30], 

7075-T6 , 3PB 

Kaufman  [34], 

7079-T6,  -T651 ,  DEN 

7075-T6,  -T651,  DEN 

Kaufman  [35], 

7178-T7651,  SEN 


c 

► 

* 

■ 

k 

l 

% 

i 

\ 

« 

« 

1 

■ 

Table  A.2  (con't) 

* 

i  Source,  material,  specimen 

Temp. 

Bo 

gi 

Class. 

J 

J 

do 

jj 

2048-T851,  (L-T),  CT 

0.75 

0.67 

1.20 

b-d 

i 

0.93 

1.04 

b-d 

1.20 

0.90 

b-d 

1.47 

0.66 

b 

2.00 

0.64 

b 

2048-T851,  (T-L),  CT 

0.75 

0.67 

1.44 

b 

0.83 

1.16 

b 

> 

1.00 

1.00 

b 

\ 

1.33 

0.91 

b 

% 

i 

1.50 

0.87 

b 

!  2048-T351,  (LT),  CT 

0.75 

0.67 

1.39 

d 

s 

1.00 

1.00 

d 

1 

1 

1.33 

0.83 

b-d 

' 

1.67 

0.95 

b-d 

, 

2.00 

0.95 

b-d 

. 

.* 

1.00 

0.50 

1.30 

d 

N* 

0.75 

1.48 

d 

1.00 

1.00 

d 

> 

1.25 

1.07 

d 

1.50 

0.95 

d 

0.50 

1.00 

1.00 

b-d 

1.50 

0.94 

b-d 

2.00 

0.91 

b-d 

> 

!  Read  and  Reed  [49], 

-460°F 

0.6 

0.25 

1.34 

b-d 

1  2219-T87,  CT 

0.47 

1.11 

b 

0.93 

1.01 

b 

Shannon  et  al.  [54], 

0.125 

2.0 

0.93 

d 

> 

2419,  DEN 

4.0 

0.87 

d 

0.25 

1.0 

1.00 

d 

2.0 

0.85 

d 

4.0 

0.77 

d 

0.50 

1.0 

1.00 

d 

1 

2.0 

0.85 

d 

7075-T7351 

0.125 

2.0 

0.96 

d 

4.0 

0.93 

d 

0.25 

1.0 

1.00 

d 

| 

2.0 

0.83 

d 

4.0 

0.62 

d 

75 

1  -  .  -  .  V-  A  V7.  ■"  '  ~ 


yvyvy 


Table  A. 2  (con't) 


Source,  material,  specimen 

Temp. 

Bo 

Oj 

Class 

GO* 

0.50 

1.0 

1.00 

d 

2.0 

0.64 

d 

Sullivan  and  Stoop  |60], 

0.61 

0.07 

0.80 

b-d 

7178-T6,  CCP 

0.10 

1.00 

d 

0.15 

0.88 

b-d 

0.21 

0.90 

b-d 

7079-T6,  CCP 

0.87 

0.04 

0.68 

d 

0.07 

0.78 

d 

0.29 

0.89 

d 

2014-T6,  CCP 

0.82 

0.05 

1.31 

d 

0.11 

1.31 

d 

0.31 

1.16 

d 

2219-T87,  CCP 

0.96 

0.03 

0.81 

d 

0.09 

0.86 

d 

7075-T6,  CCP 

0.69 

0.09 

1.89 

d 

0.44 

1  24 

b-d 

7475-T61,  CCP 

0.75 

0.14 

1.03 

d 

0.25 

1.02 

d 

1.04 

0.10 

0.47 

d 

0.18 

0.39 

d 

0.24 

0.45 

d 

1.4 

0.05 

4.77 

d 

0.14 

2.79 

d 

1.66 

0.04 

1.32 

d 

0.15 

1.18 

d 

1.15 

0.09 

1.76 

d 

0.22 

1.45 

d 

Sullivan  et  al.  [61], 

0.5 

0.06 

2.38 

d 

7075-T6,  CCP  (ELOX  precrack) 

0.13 

2.85 

d 

0.18 

3.04 

d 

0.20 

2.93 

d 

0.25 

2.58 

b-d 

0.40 

2.77 

d 

0.50 

2  14 

b-d 

7075-T6,  CCP  (fatigue  pre-cracked) 


0  5 


0  06  1.13 


d 


type  (#  tested) 

[in.] 

fication 

Johnson  [31], 

1.22 

0.04 

1.58 

b-d 

Ti-6AI-4V,  CTS 

0.32 

1.18 

b 

Jones  and  Brown  [32], 

0.25 

0.08 

1.59 

b-d 

Ti-6AI-6V-2Sn 

0.24 

1.27 

b 

0.52 

0.96 

b 

1.00 

1.00 

b 

2.04 

1.01 

b 

Poulose  et  al.  [47], 

0.75 

0.17 

1.13 

d 

Ti-6AI-4V 

0.33 

1.09 

b-d 

0.67 

1.03 

b-d 

1.00 

1.00 

b-d 

1.33 

0.86 

b-d 

Repko  et  al.  [50], 

76°F  0.075 

1.73 

1.10 

d 

DEN 

0.84 

1.08 

d 

0.33 

1.34 

d 

0.24 

1.54 

d 

0.13 

1.76 

d 

0.075 

0.24 

0.86 

b-d 

0.33 

1.67 

d 

0.84 

1.06 

d 

1.73 

0.83 

b-d 

0.225 

0.08 

1.04 

b-d 

0.11 

2.35 

b-d 

0.28 

1.41 

b-d 

0.58 

1.16 

b 

0.075 

0.24 

1.00 

d 

0.31 

1.00 

d 

0.84 

1.00 

d 

1.73 

1.00 

d 

0.225 

0.08 

1.00 

d 

0.10 

1.00 

d 

0.28 

1.00 

d 

0.58 

1.00 

d 

-1 10°F  0.075 

0.24 

0.96 

d 

0.31 

1.00 

d 

'tfl, 


Table  A.3  (con't) 

Source,  material,  specimen 

Temp.  B0 

^•B 

-320°F  0.075 

1.73 

0.24 

Shannon  et  al.  [54], 

0.125 

0.31 

0.84 

1.73 

2.0 

TF6A1-4V 

4.0 

0.25 

1.0 

0.50 

2.0 

4.0 

1.0 

Ti-8M0-8V 

0.25 

2.0 

1.0 

0.50 

2.0 

4.0 

1.0 

Srawley  and  Beachem  [57], 

80°  F  0.25 

2.0 

0.52 

TM3V-1 1Cr-3AI,CCP 

0.40 

Sullivan  and  Stoop  [60], 

0.55 

0  32 
0.24 
0.16 
0.08 

0.08 

Ti-4AI-3Mo-1V,  CCP 

0.11 

TM6V-2.5AI,  CCP 

0.645 

0.16 

0.23 

0.06 

Ti-13V-1 1 A-3AI,  CCP 

0.53 

0.18 

0.08 

Williams  [63], 

0.088 

0.12 

0.17 

2.86 

Ti-4AI-3Mo-1V,  CT 

1.48 

Table  A.4  (con't) 


G\ 

Source,  material,  specimen 

Temp. 

Bo 

>-B 

Class 

oo* 

3.80 

0.87 

d 

Kinloch  and  Gledhill  [37], 

68°F 

1.0 

0.4 

0.89 

b 

Nitrocellulose/nitro¬ 

2.0 

1.09 

b 

glycerine  propellant,  CT 

32°F 

1.0 

0.4 

1.07 

b 

2.0 

0.95 

b 

■4°F 

1.0 

0.4 

1.23 

b 

2.0 

0.85 

b 

-76°F 

1  .0 

0.4 

1.42 

b 

2.0 

0.80 

b 

Mindess  and  Nadeau  [40], 

1.0 

3.98 

0.86 

b 

Mortar,  3PB 

8.03 

0.80 

b 

12.0 

0.82 

b 

16.0 

0.84 

b 

20.0 

0.83 

b 

Concrete,  3PB 

1.0 

3.54 

1.08 

b 

8.03 

1.13 

b 

12.0 

1.10 

b 

16.0 

1.18 

b 

20.0 

1.07 

b 

Petrie  et.  al.  [45], 

0.5 

0.40 

1.13 

b 

Polycarbon.  sheet(Lexan  9030),  CT 

1.0 

1.06 

b 

Schmidt  and  Lutz  [53], 

1.97 

0.27 

0.99 

b-d 

Westerly  granite 

0.55 

1.04 

b-d 

1.04 

0.99 

b-d 

1.53 

1.00 

b-d 

2.04 

1.02 

b-d 

Smith  and  Chowdary  [56], 

0.238 

1.14 

1.05 

b 

Slip-cast  fused  silica,  3PB 

1.72 

1.23 

b 

1.92 

0.97 

b 

2.11 

1.15 

b 

Williams  and  Ewing  [64], 

0.055 

2.25 

0.87 

b-d 

Polymethylmethacrylate,  SEN 

4.50 

0.78 

b-d 

0.105 

1.19 

0.96 

b-d 

2.38 

0.83 

b-d 

0.15 

0.83 

0.96 

b-d 

1.66 

1.10 

b-d 
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APPENDIX  B.  COOLING  TIME  CALCULATIONS 


In  this  appendix  we  present  a  method  of  estimating  how  long  the  test  specimens  need  to 
sit  in  liquid  nitrogen  (LN2)  in  order  to  reach  equilibrium.  We  start  by  developing  a  simplified 

physical  model,  then  formally  state  the  problem  to  be  solved.  After  this  we  present  a  solution 
leading  to  a  means  of  obtaining  the  times.  We  close  by  discussing  the  reliability  of  these 
estimates. 

The  geometry  of  the  specimens  is  as  shown  earlier  in  Figure  3.1.  Computing  the  hold 
time  in  LN2  for  the  center  to  reach  the  same  temperature  as  the  surface  is  a  problem  of  heat 

conduction  in  a  solid.  Such  problems  have  been  extensively  handled  by  Carslaw  and  Jaeger 
[145].  Although  there  is  no  solution  for  this  exact  geometry  in  [145],  we  can  make  use  of  Carslaw 
and  Jaeger  by  setting  up  the  problem  at  hand  in  a  simplified  model.  The  simple  model  chosen  to 
represent  the  specimens  is  a  finite  length  hollow  cylinder,  Figure  B-1 . 

This  model  is  sized  so  that  the  inside  radius,  r,,  is  one  half  of  the  machined  notch  width,  as 
indicated  by  the  dashed  lines  in  Figure  B.1 .  The  outside  radius,  r0,  is  selected  so  that  the  model 

shares  the  same  length  ahead  of  the  uncracked  ligament.  The  thickness,  B,  is  taken  to  be  the 
same  as  for  the  fracture  specimen.  All  surfaces  are  awash  with  LN2,  and  therefore  we  assume  that 
they  come  to  LN2  temperature  virtually  immediately. t  The  left  half  of  Figure  B.l  is  the  region  of 

greatest  interest  in  the  specimen  and  specifically  we  are  concerned  with  the  temperature  along 
the  0  =  0  line,  the  probable  trajectory  of  the  fracture. 


t  This  assumption  was  confirmed  with  surface  thermocuple  measurements. 
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Figure  B.l .  Model  geometry  used  for  thermal  calculations 


The  general  problem  statement  for  this  model  geometry  proceeds  as  follows.  We  seek  to 
find  the  time  for  a  piece  of  1045  steel,  in  the  geometry  shown  in  Figure  B.l,  to  reach  equilibrium 
when  taken  from  room  temperature  and  immersed  in  a  bath  of  liquid  nitrogen.  In  order  to 
determine  this  equilibrium  time,  te,  we  need  to  find  the  axisymmetric  temperature  distribution,  T= 
T(r,z,t)  for  all  time  t  (>0)  throughout  R ,  where  R  is  the  region  defined  by 


R  5  {(r,8,z)  |  r  <  r0,  -n<  6<n,  -B/2  <  z  <  B/2}, 


and  r,  0,  z  are  cylindrical  polar  coordinates,  (Figure  B.l). 


The  temperature  distribution  is  to  satisfy: 


row 


BMW 


on  R  for  t>0,  where  V2  =  — +  —~  in  the  Laplacian  operator  and  K  is  the  conductivity;  and 

3r^  f  3r  dzd 


the  initial  and  boundary  conditions, 


f=Tr  att  =  0  on R 


T  =  T|  on  r  =  q,  r0  (-n  <  0  <7r,  -B/2  <  z  <  B/2,  0  <t)  and 


T  «  T|  on  z  =  ±B/2  (-rc  <  0  <7t,  rj  <  r  <  r0,  0  <t) 


The  solution  to  the  above  problem  can  be  found  in  essence  in  Carslaw  and  Jaeger  [145] 
Sections  7.10  and  8.4,  and  is: 


T  =  T|  +  (Tr  -  T,)  X¥ 


where  Xsit  V exp(-Ka^t)  ■  -^an-—  UslffOrl—  ancj 

4-d  n  Jo(onr.)  +  J0(anro) 

n=1 


4  V  (-1)n  f-K(2n+1)2rt2A  ,(2n+1 


Here  Jo(anr) is  a  Bessel  function  of  the  first  kind,  order  zero,  Yo(anr)  is  a  Bessel  function  of  the 
second  kind,  order  zero,  Uo(anr)  is  the  associated  Wronskien  (Uo(anr)  =  Jo(anr)Yo(anr0)  - 
J0(anro)Y0(anr))  and  an  are  the  eigen  values  of  the  problem,  these  eigen  values  are  obtained 
from  the  roots  of  Uo(anrj)  =  0. 


Since  the  time  we  seek,  that  to  reach  equilibrium,  is  long  (te»1),  the  series  solutions  of 
both  X  and  'f'  are  dominated  by  the  initial  terms.  This  domination  is  the  result  of  the  presence  of 
exponentials  raised  to  negative  powers  that  depend  upon  time.  In  the  4*  expression  the 
influence  of  the  terms  drops  off  quickly  since  the  power  of  the  exponent  changes  by  a  factor  of 
nine  in  the  first  step,  and  more  in  the  subsequent  ones.  Therefore  we  only  apply  the  first  term. 
However  the  behavior  of  the  X  expression  depends  on  how  an  varies,  and  that  variation  is  not 


obvious.  For  simplicity  at  this  time,  we  deal  with  the  first  term  in  the  X  series  and  comment  on  this 
choice  subsequently. 

Since  T  is  a  function  of  position,  we  must  identify  the  last  point  to  reach  the  desired 

temperature.  With  respect  to  the  z  dimension,  from  symmetry,  we  can  see  that  this  occurs  on  the 

z=0  plane.  However,  it  is  not  so  obvious  in  the  r  direction.  To  get  the  r  position  which  has  the 

highest  temperature,  we  take  the  derivative  of  (B.3),  with  respect  to  r,  and  equate  it  to  zero  to 

dT  dX  „  auo(or)  „ 

obtain  an  r ,  rm,  which  maximizes  the  temperature,  i.e. ,  we  set  — =  0  — =  0  => — - =  0. 

dr  dr  dr 

Combining  the  identity  for  Uo  and  the  expressions  3Jo(anr)/3r  =  -anJi(anr)  and  3Yo(anr)/dr  =  - 
anYi(anr)  then  yields. 

J0(«nro)Yi(anrm)  -  Jl(anrm)Yo(anr0)  =  0.  (B.4) 

Equation  (B.3)  can  be  solved  numerically  to  provide  rm.  Using  the  location  of  the  last  point  to  cool 
in  conjunction  with  truncating  the  exponential  series  in  t  then  gives  (for  t>>1),  for  the  maximum 
temperature,  Tm, 

Tm  ■  T,  ,  4(Tf  ■  T,{  ■*><■»  rj>  ,ml  1  e»p(-KleC?  ^2))  (B.5) 

LJo(ain)  +  Jo(airo)J  1  ° 


This  can  be  rearranged  to  give: 


Since  there  are  eleven  combinations  of  B's  and  r*s,  a  table  of  values  is  in  order  and  is  shown  in 
Table  B.l. 

We  note  in  Table  B.l  that  if  the  equilibrium  requirement  is  made  more  stringent  so  that  it 
requires  uniformity  to  within  say  0.1  percent,  the  cooling  times  are  increased  by  a  factor  of  about 
1 .4.  We  also  comment  upon  the  values  of  properties  used  to  compute  K,  namely  k  the  thermal 
conductivity,  p  the  density,  and  C  the  specific  heat.  These  values  are  commonly  given  in 
references  at  room  temperature  or  higher.  AJI  these  properties  vary  with  temperature,  although  k 
and  p  do  so  to  a  small  extent.  However,  C  can  change  and  usually  gets  smaller,  thus  increasing 
the  value  of  K.  As  seen  in  (B.6),  t  is  inversely  proportional  to  K,  so  that  an  increase  in  K  will 
shorten  te.  Finally  we  remark  on  the  use  of  only  the  first  term  in  the  X  series.  After  computing  the 

series  with  and  without  the  second  through  fourth  terms  for  the  two  inch  diameter  specimens,  we 
find  that  for  te=  1 00  seconds  the  value  of  the  X  expression  is  reduced  at  rm  by  about  10%. 

However  when  te  exceeds  500  seconds,  the  reduction  is  less  than  0.5%.  Hence  there  is  little 
effect.  The  additional  terms  reduce  the  value  of  X,  and  hence  te,  because  Uo(otnrn)  is  negative 

for  n=  2  and  3.  Thus  the  times  calculated  from  (B.6)  can  reasonably  be  expected  to  be  used  as 
upper  bounds. 


TABLE  B.1  Specimen  cooling  times 
Specimen  B  r0 

size  (in)  (in) 

n 

(in) 

ai 

rm 

(in) 

1  e 
(min) 

lx  3/8 

0.375 

0.421 

0.031 

7.5787 

0.182 

6.4 

lx  3/4 

0.75 

- 

- 

- 

« 

10.9 

2x  3/32 

0.093 

0.842 

0.062 

3.7894 

0.363 

0.7 

2  x3/16 

0.186 

M 

» 

- 

2.7 

2x3/8 

0.375 

M 

- 

9.6 

2x3/4 

0.75 

« 

- 

- 

" 

25.5 

2x3/2 

1.5 

- 

M 

N 

- 

43.4 

2x2  1/4 

2.25 

M 

H 

N 

M 

49.9 

4x3/8 

0.375 

1.684 

0.124 

1.8947 

0.727 

11.1 

4x3/4 

0.75 

" 

" 

" 

38.5 

4x3/2 

1.5 

- 

- 

- 

N 

102.0 

In  arriving  at  the  above  the  following  parameter  values,  which  are  independent  of  size,  are  used 
a-|r0=  3.191  arj=  0.2349  U0(airm)=  0.2833 

.Jg(ain)  Ug(airm)  l=0418Q 
Jo(ain)  +  Jo(air0) 

(T  -  T|)/(Tr-  T|)  =  0.39/390.04  °R/°R  =  0.001 
K  =  k/(pC)  =  1 .52  x  10  -4  in2/sec  (°R/°R) 


( 
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t  A  indicates  testing  after  heat  treating 


in  co  cm 

CO  CM  1^-  •*- 

coT-o>cocDcoin^r 

CM  t—  CM  CO 

Ol^tON 

■*—  CO  t-  O) 

O)  CD  (O  CM  Tf  ID  tr  T- 

CD  t-  CM  CO 

Olf  OO 

C\i  T-  T- 

COCOCMCO'M-COCO'f 

cb  cb  cb 

CM  CO  CD  CO 

CO  CO  CO  CO 

cocococococococo 

co  co  co  co 

Ui 

Ui 

LU 

ill 

> 

> 

> 

> 

< 

< 

< 

< 

If)  CO  Tf  -J 

SNNn 

CO»-0)NOOO)-i 

co  *-  co  — , 

n^fNo 

Is  O)  CD  Q 

com-c-mcnincoQ 

in  oo  q 

in  in  in  < 

CM  CM  CM  < 

m  in  m  < 

CD’tN 

co  co 

r —  r~-  r~-  r- 

oo  oo  oo 

0)0)0) 

0)0  0 

0  0  0  0  0  0)0 

•M-  ■M-  'M- 

CM  CM  CM 

Osi  CM  CM 

CM  CM  CM  CM  CM  CM  CM 

t*  ^  t* 

CM  n  O) 

Tf  m  o 

COCOCOOIOCON 

M  00  t- 

CO  CO  IO 

co  m  to 

co  co  r>-  in  co  co  co 

CO  h-  CO 

s  o  o  o 
555 

in  in  m 
r-  r~. 

o  ©  o 

n  n  n  s  f*1*  N 

CO  CO  CO  CO  CO  CO  CO 

ddddddo 

in  in  in 

CM  CM  CM 

CM  CM  CM 

42  1 

4  2  2 

423 

44  1 

442 

443 

«-  CM  (O  f  If)  CO 

OO  OO  00  00  OO  00  00 

■M-  ^  ^ 

2  1.3  1 

2  1.3  2 

2  1.3  3 

(con't) 


I 


TABLE  C.3  Details  of  cyclic  experiments. 


